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XENON10 is the first two-phase xenon time projection chamber (TPC) developed within the 
XENON dark matter search program. The TPC, with an active liquid xenon (LXe) mass of about 
14 kg, was installed at the Gran Sasso underground laboratory (LNGS) in Italy, and operated 
for more than one year, with excellent stability and performance. Results from a dark matter 
search with XENON10 have been published elsewhere. In this paper, we summarize the design 
and performance of the detector and its subsystems, based on calibration data using sources of 
gamma-rays and neutrons as well as background and Monte Carlo simulations data. The results 
on the detector's energy threshold, energy and position resolution, and overall efficiency show a 
performance that exceeds design specifications, in view of the very low energy threshold achieved 
(<10 keVr) and the excellent energy resolution achieved by combining the ionization and scintillation 
signals, detected simultaneously. 

PACS numbers: 95.35. +d, 29.40.Mc, 95.55.Vj 
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I. INTRODUCTION 

The XENON10 detector, a two-phase time projection 
chamber (TPC) containing 25 kg of pure liquid xenon 
(LXe), was built and operated as the first practical pro- 
totype of the XENON dark matter search program pQ. 
The R&D phase which culminated with the XENON10 
TPC, involved several smaller size detectors, dedicated to 
study many of the properties and performance character- 
istics of a two-phase TPC in the energy regime relevant 
to the direct detection of dark matter particles scattering 
off Xe nuclei. The most important results obtained with 
these prototypes are published [2H5]. 

In addition, several technical developments carried out 
during this phase played an important role in establish- 
ing the feasibility of the XENON approach for dark mat- 
ter detection. Among these, two are particularly impor- 
tant: a) the development by Hamamatsu Photonics of 
compact metal channel photomultipliers (PMTs) for the 
detection of the LXe scintillation light, with continuous 
improvement in quantum efficiency and radio-purity [6]- 
[8] ; b) the development of a pulse tube refrigerator (PTR) 
optimized for LXe temperature [9 to achieve the required 
long-term stability of a cryogenics detector filled with a 
large LXe volume. 

The goal of the XENON 10 experiment was to demon- 
strate the achievable energy threshold, background rejec- 
tion power and operation stability of a two-phase TPC at 
the 10 kg fiducial mass scale, prior to the realization of a 
100 kg scale detector, originally proposed as unit module 
for a ton scale XENON dark matter search. 

The fast-paced R&D effort did not emphasize mate- 
rials radio-purity and XENON 10 was built with largely 
off-the-shelf components, not screened or selected to min- 
imize backgrounds. The TPC and its associated cryo- 
genics, purification and data acquisition (DAQ) systems 
were developed and tested at the Nevis Laboratories of 
Columbia University prior to installation at the Ital- 
ian Gran Sasso Underground Laboratory (LNGS) [ID] . 
The commissioning of the XENON10 detector at LNGS 
started in April 2006. Following the completion of the 
shield system, source calibration and background data 
were acquired for several months, in stable conditions. 
A blind analysis of 58.4 live days of data, acquired be- 
tween October 6, 2006 and February 14, 2007, and us- 
ing a fiducial mass of 5.4 kg, excluded previously unex- 
plored parameter space, setting improved upper limits 
on both spin-independent [IT] and spin-dependent [12] 
cross- section for dark matter scattering off nucleons. 

The XENON10 results, and the unusually fast time 
scale in which they were achieved, validated the scien- 
tific reach of a position sensitive, homogeneous and self- 
triggered LXeTPC for dark matter direct detection. New 
and improved LXe dark matter detectors, with a target 
mass at the 100 kg scale, are currently in operation [15] 



or under construction [16] [17] . They promise to advance 
the field with significant improvement in sensitivity over 
the next few years. In this paper we wish to summa- 
rize the development and performance of the XENON 10 
detector and associated systems. The plan of the pa- 
per is the following. In Section 2 we describe the design 
of the detector, cryogenics, purification and slow-control 
systems, followed by a description of the shield and data 
acquisition systems. In Section 3 we discuss the data 
processing, the algorithms and cuts used to infer the pa- 
rameters relevant for a dark matter search. In Section 4, 
we present results from gamma and neutron calibration 
of the XENON10 TPC. Finally in Section 5, details on 
the analysis of XENON 10 data used for the published 
dark matter search results are discussed. The appendix 
to the paper deals with a summary of the sources of 
backgrounds in the XENON 10 experiment, in compar- 
ison with Monte Carlo predictions. 

II. THE XENON10 DETECTOR 
A. Principle of Operation 

The principle of operation of the XENON 10 two-phase 
(liquid-gas) time projection chamber (TPC) is shown 
schematically in Fig. [I] The energy deposited by a parti- 
cle interaction within the active LXe volume is detected 
by the simultaneous measurement of ionization electrons 
and of primary scintillation photons (£1), with a wave- 
length of 178 nm [13], produced by the de-excitation 
to the ground state of excited diatomic Xe molecules 
(XeJ) [H]. Both direct excitation of atoms and electron- 
ion recombination lead to the formation of XeJ . 

Under the application of an electric field, Ed, on the or- 
der of 1 kV/cm, the electrons which escape recombination 
with the parent positive ions, drift towards the liquid-gas 
interface with a velocity of about 2 mm//is. Once they 
reach the liquid surface, they are extracted into the gas, 
where a much stronger electric field, £ e , accelerate them 
leading to secondary scintillation photons emission (S2). 
In the XENON10 TPC, photomultiplier tubes (PMTs), 
coupled directly to the sensitive liquid and gas regions, 
are used to detect both SI and S2 signals. As shown in 
[2], the extraction yield is 100% for a field greater than 
10 kV/cm. The number of secondary scintillation pho- 
tons emitted in the gas is proportional to the number of 
electrons, and thus to the energy deposited in the liq- 
uid by the incoming particle. We refer to this signal as 
proportional scintillation signal. The following equation 
expresses the dependence of this signal on the operation 
conditions: 

N 7 = aN e (e e /p- /3)px. (1) 

Here p is the gas pressure and x is the distance traveled 
by the electrons in the gas which largely determines the 
width of the S2 pulse. The units for p, e e and x are 
bar, kV/cm and cm, respectively. N e is the number of 
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electrons extracted from the liquid to the gas phase, iV 7 
is the number of proportional scintillation photons pro- 
duced, a (the amplification factor) and f3 (the threshold 
of the reduced field for proportional light production) 
are experimentally determined values. For the case of 
Xe, the reported values for a vary from 70 to 140 and 
p values are between 0.8 and 1.0 kV/cm/bar [2| [18} fT9] . 
The number of secondary photons generated by one drift- 
ing electron is sufficiently large to be detected by PMTs, 
thus a two-phase detector is single-electron sensitive. In 
XENON10 operating conditions, the values of p, e e and x 
were 2.1 bar, 12 kV/cm and 0.25 cm, respectively, yield- 
ing ~100 photons per extracted electron. 

The different amplitude of the SI and S2 signals asso- 
ciated with nuclear recoils, such as produced by WIMP 
or neutron interactions, and with electron recoils, such as 
produced by background beta and gamma-rays, provides 
the basis for background discrimination in a two phase 
TPC such as XENON10. Since electron diffusion in LXe 
is small [20], the proportional scintillation pulse is pro- 
duced in a small volume with the same XY coordinates 
as the interaction site, allowing 2D localization. The Z 
coordinate is inferred from the drift time measurement, 
and the known electron drift velocity at the operating 
field. The TPC approach and the strong self-shielding 
property of LXe (p w 3 g/cc, Z = 54), make the XENON 
concept very powerful for background identification and 
rejection based on event topology: multiple scatter events 
from the overwhelming majority of gamma rays and from 
neutrons can be identified and rejected. The sensitive 
target can be "fiducialized" to keep only the inner core 
free of background. 

The capability to accurately identify the 3D position 
on an event-by-event basis, also allows to correct the po- 
sition dependence of both the direct and proportional 
scintillation signals, which in turn results in improved 
energy resolution. For a recent review of the properties 
of LXe and its response to radiation we refer to [20] . 

One of the challenges of a two-phase detector is to ef- 
ficiently detect the small number of primary scintillation 
photons associated with the low energy events of interest 
in a dark matter search. Because of the large refractive 
index of LXe 2T] and consequent total internal reflection 
at the liquid/gas interface, the presence of PMTs at the 
bottom of the active drift volume is essential for efficient 
SI light collection. 

When the XENON concept was proposed, the devel- 
opment of PMTs capable to operate directly in LXe had 
just began. The typical quantum efficiency of these early 
PMTs at 178 nm was in the 10% range. To achieve 
high light collection of the 51 signal and thus a low 
energy threshold required for a dark matter search, the 
baseline XENON detector design used a Csl photocath- 
ode in place of a common cathode, and photomultipliers 
or other type of photodetectors in the gas region [22] . 
During the XENON R&D phase, different photodetec- 
tors were tested [23 , 24 and also charge-sensitive readouts 
based on the Gas Electron Multiplier concept [25 . The 





~1 \is width 




Anode Grid 



Cathode Grid 



FIG. 1: (Color online) Principle of operation of the XENON10 
two-phase TPC. A particle interaction in the liquid produces 
primary scintillation light (SI) and ionization electrons. The 
electrons drift under an electric field Sd (about 1 kV/cm) un- 
til they are extracted from the liquid into the gas where they 
are accelerated by an electric field s e (about 10 kV/cm) pro- 
ducing proportional scintillation light (S2). Two arrays of 
photomultipliers, one in the liquid and one in the gas, detect 
simultaneously the SI and S2 light signals. 



operation of Csl photocathodes to detect the direct light 
in two phase prototypes confirmed [26 the high QE first 
measured by the Columbia group [27]. However, the com- 
plexity of operating a TPC under the very strong electric 
field required by the Csl photocathode, while suppress- 
ing positive photon feedback, and the fast improvement 
of the performance of Hamamatsu PMTs for LXe opera- 
tion, led us to abandon the Csl photocathode approach. 
Prior to adopting the Hamamatsu one-inch square R8520 
PMT for XENON10, a two-inch diameter R9288 PMT 
was extensively tested [2j [26] . 



B. Electrodes Assembly and Electric Fields 

Fig. [2] shows a 2D mechanical drawing of the 
XENON10 detector assembly, and Fig. [4] shows a pic- 
ture of the detector mounted on the movable door of the 



passive shield (see section III). A closer view of the in- 
ner TPC structure is shown in Fig. [3] The TPC active 
volume is defined by a polytetrafluoroethylene (PTFE) 
cylinder with an inner diameter of 20 cm and a height 
of 15 cm for a total active mass of about 14 kg of LXe. 
PTFE is used as UV light reflector [28] and as electrically 
insulating support for the TPC structure. 

To produce the electric fields in the liquid and gas re- 
gions, four wire meshes, also referred to as grids, are used; 
two in the liquid (cathode and gate grid) and two in the 
gas (anode and top grid). The 0.203 mm thick electro- 
formed meshes, are made of electropolished 304 stain- 
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less steel,with a bar width of 0.182 mm and 2.0 mm x 
2.0 mm square holes, insuring good optical transmission. 
The cathode, at the bottom of the PTFE cylinder, and 
the gate grid define the TPC 15 cm drift distance. Field 
shaping rings, made of 0.5 mm thick copper and spaced 
by 0.76 cm, are mounted outside the PTFE cylinder to 
insure a uniform electric field across the drift volume. 

A 5 mm gap separates the gate grid from the anode 
and the anode from the top grid. The liquid level is 
between anode and gate grid and determines the extrac- 
tion field. The liquid level, 2.5 mm above the gate grid, is 
kept constant by the use of a pressurized cylinder closed 
on the top, similar to a diving bell. The pressure is pro- 
vided by the return of th e gas circulating through the 
purifier (see Section HE). Custom-made capacitors are 
used to monitor the liquid level: one cylindrical capaci- 
tor is used to measure the liquid level and four parallel- 
plate capacitors are used to measure the inclination of 
the detector, to ensure that the liquid level is parallel to 
the grids and thus that the extraction field is uniform. 
One of the parallel-plate capacitors is filled with PTFE 
and is used as a reference capacitor: the capacitances 
of other parallel-plate capacitors are compared with the 
reference value to control the inclination of the detector. 
The XENON10 cryostat was equipped with leveling feet 
which could be adjusted from outside the shield structure 
to achieve the required degree of leveling. 

To minimize passive LXe outside the TPC, the PTFE 
cylinder is surrounded by an outer PTFE cylinder, with 
cut-outs for the resistors of the voltage divider network, 
directly mounted on the field shaping wires. A photo of 
the assembled TPC structure, with the top three meshes 
clearly visible, is shown in Figure [5] For most of the 
data acquired with XENON10, the top and gate grids 
were biased to -1.15 kV, with the anode and cathode at 
+3.5 kV and -13 kV, respectively. With these voltages, 
the drift field in the liquid was 0.73 kV/cm while the 
field in the gas was 12 kV/cm. The high voltage for the 
gate, anode, and top meshes was provided by CAEN 1733 
and 1833 power supplies [29 , while the higher voltage 
for the cathode was provided by a Heinzinger power sup- 
ply [30] . Commercial HV feedthroughs and cables carried 
these voltages to the cryostat. To carry the high voltages 
from the vacuum of the cryostat to the cathode mesh, a 
custom-made PTFE insulated feedthrough was used. In- 
side the TPC structure and below the LXe level, bare 
wires were used to carry the HV to the upper meshes, 
taking advantage of the excellent dielectric properties of 
LXe and of PTFE. 

Results from electric fields simulations for the 
XENON10 TPC are shown in Fig. § for a drift field 
of 1 kV/cm. The field is uniform up to 3 mm from the 
PTFE wall. For the dark matter search results, only 
events with at least 10 mm distance from the edge of 
the active volume were accepted. Gamma-ray calibra- 
tion data confirmed that no charge was lost for events in 
this region. On the bottom of the detector, a region of re- 
versed field direction exists between the cathode and the 
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FIG. 2: (Color online) Schematic drawing of the XENON10 
detector. 
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FIG. 3: (Color online) A close-up view of the XENON10 TPC 
structure. 



bottom PMT array, separated by a distance of 1.2 cm. 
The field reversal extends slightly into the sensitive vol- 
ume, affecting a region less than 0.9 mm above the cath- 
ode. 

The same set of field simulations was used to estab- 
lish that the displacement of the true XY event position 
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FIG. 4: (Color online) Photo of the XENON10 detector and 
shield, with the shield door open. 




FIG. 5: (Color online) Photo of the assembled XENON10 
TPC structure. 



and the mean XY position from the proportional scin- 
tillation was negligible (< 0.6 mm). The impact of the 
field distribution in the gas on the energy resolution of 
the detector was also studied. For the inner region of the 
meshes, it was found that the relative path length vari- 
ation of electrons traversing the gap from the liquid-gas 
surface to the anode mesh is < 20% FWHM. 



C. PMTs and Calibration System 

The active LXe volume is viewed by 89 PMTs: the 
bottom array of 41 PMTs is located 1.2 cm below the 
cathode, fully immersed in LXe. The 48 PMTs of the 
top array are in the gas. The XENON 10 PMTs are 1" 
x 1" square Hamamatsu R8520-06-AL [31] designed to 
work in LXe and optimized for low radioactivity. They 
have a bialkali photocathode and a quartz window, with 
a typical quantum efficiency > 20% at 178 nm. They are 




FIG. 6: (Color online) Simulated electrostatic field in 
XENON10. Arrows and colors indicate the field direction 
and strength, respectively. The field cage was optimized for 
field uniformity in the sensitive region. The region below the 
cathode shows a strong reversed field. 



compact (only 3.5 cm tall), metal channel PMTs, with 
10 multiplication stages and a total amplification of a 
few 10 6 . They are certified for operation at -100° C and 
up to 5 bar. The radioactivity of these PMTs was mea- 
sured with a ultra sensitive HPGe detector [32] at the 
Gran Sasso Low Background facility [33]. Their contam- 
ination in 238 U/ 232 Th/ 40 K/ 60 Co was measured as low as 
(0.25 ± 0.04)/(0.21 ± 0/05)/(9.3 ± l.l)/(0.59 ± 0.05) 
mBq/PMT. 

A photo of the top array, mounted inside the "diving 
bell" structure which defines the liquid level, is shown 
in Fig. [7] The custom-made cathode HV feedthrough 
and the level meter are also visible. The high voltage 
for the PMTs is provided by CAEN 1833 power supply 
boards [29]; this system featuring single channel control, 
allowed to easily equalize the gain of the PMTs, by ad- 
justing individual voltages. 

The HV divider for the PMTs was mounted on a Cirlex 
substrate, with surface mount components. The Cirlex 
boards for each of the two arrays are mounted on cop- 
per plates with PTFE spacers for insulation. Despite the 
small mass, the ceramic of the resistors and capacitors 
contribute the largest fraction of the radioactive back- 
ground of the HV divider boards. Heating of the LXe 
by the resistive dividers is kept to a minimum by choos- 
ing high resistance values (10 MSI), allowed because of 
the low event rate and low light level during dark matter 
search data-taking. The typical power consumption of 
the XENON10 PMTs was about 5 mW per PMT. 
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Negative HV was applied to the cathode of each PMT, 
via unshielded Kapton wire. The average HV was about 
750 V, with a maximum of 900 V. The signal was carried 
out with 50 Q teflon cables (RG178U), stripped of their 
shielding, to avoid trapped air which would impact the 
liquid purity. PMTs signal and HV cables were carried 
outside by 48 pin Burndy feedthroughs [3l] . The radioac- 
tivity of these feedthroughs was reduced by replacing the 
ceramic plugs holding the pins with plugs made out of 
polyether ether ketone (PEEK). 




FIG. 7: (Color online) XENON10 top PMT array, with the 
cathode HV feedthrough and liquid level meter also visible. 

The 89 PMTs were routinely calibrated with the goal 
of equalizing the gain and to correct for differences in 
the quantum efficiency. The gain equalization was per- 
formed using light emitting diodes (LEDs). Two LEDs 
were used: one was mounted close to the top array, to 
equalize the bottom PMTs; the other was mounted close 
to the bottom array, to equalize the top PMTs. Both 
LEDs were covered with a PTFE cap acting as a dif- 
fuser to ensure a uniform illumination of all the PMTs in 
each array. The LEDs were driven by a pulse generator, 
which synchronously activated also the main DAQ trig- 
ger: the pulse duration was 6 /is, but the charge on each 
channel was integrated only in a 1 /is window, 4 /as after 
the activation of the LED, to minimize electronics noise 
generated when the LED was switched on. 

The gain measurement was performed in single photo- 
electron (p.e.) regime. Under the assumption that the 
p.e. spectrum is the sum of a Gaussian noise peak and a 
Gaussian single p.e. peak, each with a Poisson-like dis- 
tribution, and with a single photoelectron to noise ratio 
of 1/5, less than 2% of the events had 2 or more pho- 
toelectrons (a typical single p.e. spectrum is shown in 
figure [§. To ensure that the single p.e. condition was 
met for all the PMTs, the calibration data were acquired 
under different illumination conditions. The spectra for 
each PMT were fitted with a function which included 
the noise Gaussian, the single p.e. Gaussian and a third 
Gaussian accounting for the multiple p.e. contribution. 
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FIG. 8: Single photoelectron spectrum of a XENON10 top 
PMT, acquired during a routine LED calibration. Clearly 
visible is also the noise peak. 



After measuring the gain for each channel, the ampli- 
fication factor was equalized for all PMTs, by adjusting 
the single PMT operating voltage. After the equaliza- 
tion, the gains were measured again and the resulting 
amplification factors were then fed to the XENON10 data 
analysis programs (see Sec. [V]) and used as a conversion 
factor between the measured charge and the actual p.e. 
number. During normal operation, such a gain measure- 
ment was repeated at least once a week, to ensure the 
stability of the PMTs gain. The time evolution of the 
gain, measured during the WIMP search data taking, is 
shown in Fig. [9] for two typical PMTs. The variation in 
PMT gain was < 2%. 
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FIG. 9: Stability of the gain of two typical XENON10 PMTs 
during the WIMP search data taking period. 



D. Cryogenic System 

A reliable and stable cryogenics system was an es- 
sential requirement for the XENON 10 experiment since 
both PMTs gain and the proportional light yield vary 
with temperature. The XENON10 TPC vessel contain- 
ing the 25 kg of LXe, made of stainless steel (SS), was 
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wrapped with aluminized mylar and surrounded by a vac- 
uum cryostat, also made of SS, with typical vacuum in 
the mTorr range. The cryogenics system was based on 
a Pulse Tube Refrigerator (PTR) with 100 W of cooling 
power at 165 K, with a 3.5 KW compressor. This type of 
cryocooler was developed specifically for Xe liquefaction 
and re-condensation in the LXe scintillating calorimeter 
of the MEG experiment [9]. To meet the ultra-high pu- 
rity requirement for the XENON 10 TPC, the cryocooler 
was placed outside of the detector vessel, mounted on 
the vacuum cryostat top flange and attached to a copper 
plate which worked as re-condenser for the evaporated Xe 
gas. The cooling power was sufficient to cool-down the 
detector, liquefy the gas and re-condense it. The desired 
temperature of the LXe was adjusted by a Proportional- 
Integral-Derivative (PID) controlled heating element on 
the cold head of the PTR. Due to the narrow temper- 
ature margin of less than 4 degrees between the liquid 
and solid phase of Xe, temperature control during liq- 
uefaction is especially important. The typical operating 
temperature for XENON10 was -94° C with a stability 
better than 0.05° C. At this temperature the Xe vapor 
pressure is ~2.1 atm. With this cryogenic system we 
have been able to achieve the level of long term stability 
required for a dark matter experiment. Fig. [To] shows 
the stability of the LXe temperature and corresponding 
Xe vapor pressure during several months of continuous 
operation of XENON10. 
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FIG. 10: Stability of the liquid Xe temperature and corre- 
sponding vapor pressure in the XENON10 detector during the 
calibration and WIMP search data-taking period. 

The XENON 10 cryogenics system included an inde- 
pendent LN2 system for emergency cooling. A pressure 
rise in the detector, above a set value, would automat- 
ically trigger a flow of LN 2 through a cooling coil next 
to the PTR. The TPC pressure sensor and the solenoid 
valve control were powered by a dedicated uninterrupt- 
ible power supply. The emergency cooling system would 
keep detector in a safe sate even in the extreme case of 
power failure and no access to the laboratory for more 
than 24 hours. We note that the large Xe mass and high 
thermal capacity stabilizes the temperature and the pres- 



sure against abrupt changes. A rupture disk would finally 
burst if the absolute pressure in the detector would rise 
above 3.5 atm, to avoid damage to the PMTs. 



E. Xe Purification and Recirculation System 

The 15 cm electron drift gap in the XENON10 detec- 
tor imposes ultra-high purity requirements on the LXe, in 
order to minimize charge loss by electron attachment to 
impurities molecules. The overall concentration of elec- 
tronegative impurities (Oxygen equivalent) must be well 
below 1 part per billion (ppb). Moreover, stable oper- 
ation of the detector demands that this high purity be 
maintained over long periods of time. For detector and 
gas system ultra high vacuum materials and procedures 
were used, and for cleaning the Xe gas a single purifier, 
a high temperature SAES getter (PS4-MT3-R) [35] was 
used. 

However, the PMTs limited the bake-out temperature 
of the detector to 70° C. In addition, LXe is a good sol- 
vent due to Van der Waals interactions with impurity 
molecules. Thus, once pure Xe is liquefied in the detec- 
tor, the purity level can easily deteriorate. We therefore 
implemented a closed circulation system to continuously 
clean the gas through the getter, after initial liquefaction. 
Such a system was successfully developed and tested at 
Columbia on several smaller scale prototypes, prior to 
adapting it to XENON10. For gas circulation a double- 
diaphragm pump (KNF-N143.12E) [36] was used. With 
a gas flow rate of 2.6 slpm, limited largely by the PTR 
cooling power, it would take about one month to reach 
the purity level required for negligible electron attenu- 
ation. A schematic of the gas purification and recircu- 
lation system used on XENON10 at LNGS is shown in 



Fig. 11 
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FIG. 11: Schematic of the Xe purification system with closed- 
loop re-circulation. 



8 



F. Slow Control System 

A Slow Control System (SCS) was developed to mon- 
itor all essential run-time parameters of the XENON10 
experiment, controlling the status of various hardware 
components, triggering alarms and transmitting impor- 
tant parameters. 

The SCS monitors over 290 parameters, which include 
gas pressure, cryogenic temperatures, flow rate, liquid 
level, grids and PMTs high voltages, DAQ trigger and 
acquisition rate, room temperatures and detector inclina- 
tion. Figures [To] is an example of the detector parameters 
monitored by the SCS during the dark matter run. 

The SCS consists of 4 parts: (i) server, (ii) monitor 
client, (iii) alarm client, and (iv) history plotter. The 
server establishes communication with all the different in- 
struments monitored, makes all the parameters available 
for the clients, and stores the information to disk. The 
monitor client allows each user access to all the parame- 
ters monitored over the last 24 hours. The alarm client 
triggers the alarm system in case one of the preassigned 
parameters falls outside the allowed range. Alarms are 
sent via email and as text messages to the cell phones of 
the appropriate personnel. In addition to this automated 
system, the LNGS personnel continuously monitors the 
gas pressure. The history plotter allows the user to ac- 
cess the information stored for each parameter for any 
particular time. 

The SCS is a platform independent software devel- 
oped exclusively for the XENON 10 experiment using the 
Java programming language. The communication be- 
tween the different SCS components is done with Java 
Remote Method Invocation (RMI). 



III. THE XENON10 SHIELD 

The XENON 10 detector is protected from external 
background by a cubic steel-framed structure, consist- 
ing of 20 cm high-density polyethylene (HDPE) inside of 
20 cm Pb; a schematic is shown in Fig. [12] The Pb was 
supplied in 5 x 10 x 20 cm bricks and was stacked so as 
to avoid any line-of-sight penetration along the cracks. 
The outer 15 cm of Pb have an activity in 210 Pb of 
about 560 ± 90 Bq/kg, while the inner 5 cm are a low- 
background Pb obtained from Fonderies de Gentilly [37] . 
It was measured at a germanium counting facility [32j [33] 
to have an activity in 210 Pb of (17 ± 5) Bq/kg. Specific 
activities for various shield components are shown in Ta- 
ble m 

A Monte Carlo simulation shows that 20 cm Pb results 
in an attenuation of the external gamma flux larger than 
10 5 while its internal activity leads to a contribution of 
less than 0.32 cts/keVee/kg/day (dru) (E < 25 keVee) 
to the raw event rate inside the shield cavity. With the 
additional self-shielding of the outer 2 cm of Xe (as used 
in the WIMP search data analysis [11), the contribution 
of the Pb activity to the electron recoil background drops 




FIG. 12: (Color online) The XENON10 shield structure. 
Common and low-radioactivity Pb bricks are shown in yel- 
low and purple respectively, while polyethylene is light blue. 
Not shown are the 20cm of polyethylene on the "door" and 
the 15cm below the shield structure. 



below 0.05 dru (E < 25 keVee), becoming sub-dominant 
(see Fig. [13). 



Contribution of Pb-21 and Rn-222 to the XENON1 background 
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FIG. 13: (Color online) Simulated contributions of 210 Pb (17 
Bq/kg) and 222 Rn (5 Bq/m 3 ) to the XENON10 background, 
for the entire 14 kg active volume (red, blue) and also for a 
5.6 kg fiducial volume (cyan, magenta); the fiducial volume 
corresponds to 2 cm Xe self-shielding (radial) and 3 cm each 
(Top/Bottom). 

The 20 cm of HDPE provides reduction of external 
fast-neutron flux by a factor of 90. Placing the Pb outside 
the HDPE has the benefit that the neutron scattering on 
the Pb reduces the energy of the neutrons, thus increas- 
ing the efficacy of the HDPE. The dominant source of 
neutrons is expected from (a, n) and fission reactions in 
the surrounding rock, as well as from cosmogenic produc- 
tion in the rock and Pb. The ambient neutron flux from 
natural radioactivity in the rock has been measured to 
be about 4 x 10" 6 n/cm 2 /s [38j|39], and from this the in- 
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tegrated rate of neutron interactions in XENON 10 is cal- 
culated to be < 0.015 neutron/month in the dark matter 
region of interest. The muon flux in the cavern was mea- 
sured to be 1.06 ± 0.03 muon /m 2 /hour [40]; from this, 
and the muon-induced neutron yield in rock and in Pb, 
the expected rate for cosmogenic neutrons is expected to 
be about 13 times higher [41] than from natural radioac- 
tivity. 

The internal shield cavity is 107.5x90x90 cm, of which 
the detector occupies about 30% of the volume. Care 
was taken not to leave any trapped air pockets between 
the lead and polyethylene, or anywhere else in the struc- 
ture. An additional 15 cm HDPE is also present be- 
low the shield structure. All seams were coated with a 
low-radioactivity silicon-gel, and a 10 mm thick rubber 
seal runs in a 5 mm deep trench along the perimeter 
of the polyethylene on the door; thus the cavity is able 
to sustain a slight over-pressure, and boil-off N2 from a 
low-pressure dewar is passed continuously into the top 
of the cavity at a rate of 1.5 ± 0.1 standard liter /min, 
and allowed to exit at the bottom. The concentration of 
222 Rn in the shield cavity, typically about 130 Bq/m 3 , 
was reduced by the N2 purge to a value of < 5.5 Bq/m 3 
(consistent with the sensitivity floor of the RAD- 7 radon 
measuring device [42 J. The purge process occurred in 
about 30 hours. The expected contribution to the elec- 
tron recoil background (E < 25 keVee) from this level of 
Rn activity is 7 x 10 -3 dru in the inner fiducial volume 
(Fig. |13| . 

The shield structure is secured by steel panels along 
the outer walls, with bolts passing through the HDPE 
sheets. The XENON10 detector is supported by stainless 
steel I-beams on a single wall, which also serves as a door 
and can be moved laterally on rails. All cryogenic and 
electronic feed-throughs pass through this wall, using a 
Z-shape where possible, and as far off-center as possible, 
to minimize line-of-sight. The overall diameter of each 
opening was also minimized, and gaps in the Pb (after 
placing the cables) were filled with low-radioactivity Pb 
shavings. Silicon-gel was used to fill holes in the HDPE, 
also ensuring air-tightness in these locations. During the 
construction all components were cleaned and prepped 
with Ethanol before mounting them on the shield struc- 
ture, no cables or pipes have to be disconnected when 
opening the door to access the cryostat. 



IV. THE XENON10 DATA ACQUISITION 
SYSTEM 



To maximize the information available from time struc- 
ture and amplitude of the primary (SI) and secondary 
(S2) pulses, the signal from 89 PMTs are individu- 
ally digitized by an ADC array at 105 MHz (about 
10 ns/sample). The SI pulses from the direct scintilla- 
tion light of a typical background event in the energy re- 
gion of interest have 4-27 p.e. distributed over a width of 
<200 ns (with a characteristic decay time of ~27 ns [43 J. 



The S2 pulses from the secondary proportional scintilla- 
tion have x 200 the number of photoelectrons distributed 
over about 1 /is. The distance between the SI pulse and 
the S2 pulse vary from to 80 /is, the maximum drift 
time for electrons extracted in the liquid. Events with 
multiple scatters in the active Xe volume will have more 
than one S2 pulse. To guarantee that the SI pulse and 
all S2 pulses associated with an event are digitized, a full 
waveform of 160 /is is recorded for each event. 



30 
25 
0) 20 

Q- 

E 

CO 

is 15 

Cl 

10 




5 phe 



«- 0.20 \xs - 



2573 phe 



10 20 



30 



40 



50 60 



70 



80 90 100 



FIG. 14: (Color online) SI and S2 pulses for 2.3 keVee gamma 
event in XENON 10 background data. Signals from each PMT 
hit are plotted separately. An enlarged view of the SI pulse 
(5 p.e.) is also shown. See [10] for more details. 

The digitizers are Struck 3301 fast ADC (VME bus) 
modules [44 , with sampling frequency 105 MHz and 14- 
bit ADC resolution (+0.1 V ^ 4.9 V). The signal from 
each PMT is amplified using Phillips PS776 xlO am- 
plifiers, and a single-pole RC filter (30 MHz) at the in- 
put to the ADC prevents aliasing. The filter also has 
a slight shaping effect; because a single photo-electron 
pulse width is similar to the inverse sample rate of the 
digitizer, this improves the single photoelectron sensitiv- 
ity (see Fig. 14). An on-line software baseline-flattening 
algorithm sets small fluctuations in the baseline (< ±8 
digitizer bins) to zero, allowing for very efficient lossless 
data compression by a factor of ~20 using open source 
software (gzip). The data is not further modified or pro- 
cessed before being recorded to disk. 

The sum of the 30 central PMTs of the top array pro- 
vides the global S2-sensitive event trigger (see Fig. [l5| ). 
The summed signal is amplified x 10 by a CAEN N968 
with integration time constant 1 /is; this value was se- 
lected based on the width of the S2 signal. The shaped 
signal is then passed to a CAEN voltage discriminator. 
The S2 trigger was verified to have 100% efficiency for a 
square pulse voltage input of width 1 /is and a step size 
16 mV. Since a typical photoelectron from a PMT with 
2 x 10 6 gain has an integrated area of 16 mV-ns, this is 
equivalent to 100 photoelectrons. Considering that about 
55% of S2 scintillation is recorded by the top PMT ar- 
ray, it is inferred that the trigger efficiency is 100% for 
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TABLE I: Results from radioactive measurements for relevant shield components (values in mBq/kg). 





238 TT 




210pi 


235tt 
U 


40 K 


137 Cs 




54 Mn 


Pb (inner) 


< 3.9 


< 6.8 


(17 ±5) xlO* 


< 20 


< 28 


< 0.85 


< 0.19 




Pb (outer) 


< 1.6 


< 5.7 


(560 ±90) xlO 3 


< 51 


14 ±6 


< 2.1 


< 1.1 




Polyethylene 


< 5.2 


< 6.6 




< 3.2 


< 64 


< 2.6 






Steel (I-beams) 


7.8 ± 3.2 


< 4.1 




< 3.1 


< 15 


< 1.4 


170 ±3 


3.3 ±0.6 


Steel (side panels) 


< 3.5 


< 4.7 
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< 8.3 
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< 42 


< 2.8 
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FIG. 15: (Color online) Schematics for the XENON10 Data 
Acquisition System electronics. 



S2 events of 182 phe or greater. Because a single ionized 
electron extracted from LXe typically yields 24 ± 5 pho- 
toelectrons (after taking into account light collection and 
PMTs quantum efficiency), 182 photoelectrons in the S2 
signal corresponds to 7 — 8 ionized electrons. This is 
the expected charge from an event with about 1 keV nu- 
clear recoil equivalent energy [4 . The trigger signal is 
distributed from the discriminator to all ADC modules 
simultaneously using a Fan-Out chain. The S2 trigger 
provides a lower threshold than is possible with an Si- 
based trigger (as is evident from Fig. [l4| ); moreover, an 
S2-trigger allows for elimination of Sl-only events due to 
scatters in the charge-insensitive regions of the detector, 
thus reducing the trigger rate. A trigger hold-off is used 
to prevent the system from triggering on after-pulsing 
following a large S2 pulse. The hold-off requires that no 
S2 trigger candidates occur for the previous 200 /is be- 
fore allowing a trigger to start acquisition. The hold-off 
time is included in computing the acquistion dead time. 
For gamma calibration data runs, a high energy SI veto is 
added to reduce the trigger rate. The veto is based on the 
sum of the 9 center bottom array PMTs and eliminates 
events with an SI > 150 keVee. Using the S2-trigger and 
no high energy veto, the acquisition rate during WIMP 
search mode is ~2.5 Hz, leading to a dead time of 7%. 



V. THE XENON10 DATA PROCESSING AND 
ANALYSIS 

A. SI & S2 Pulse Identification 

The raw data for an event consist of 88 waveforms of 
16350 samples and is reduced to physical parameters us- 
ing two fully separate and parallel analysis chains. One 
is written in ROOT [45], and the other uses the commer- 
cial MATLAB software; both were specifically developed 
for the XENON10 data analysis. The high-level struc- 
ture of the data reduction is similar in both cases, and 
proceeds in three stages: i) preprocessing the waveforms, 
ii) searching for pulses, and finally iii) computing the re- 
duced quantities associated with each scintillation pulse. 

In the preprocessing stage, the samples of each PMT 
waveform are zeroed if they are below a threshold of ~1 
mV, which is much lower than the single photoelectron 
level and which depends on the noise characteristics of 
the channel. The waveforms of all channels are then 
added into a total waveform that is used to search for 
SI and S2 pulses. 

The pulse searching stage operates in two steps: it 
first looks for S2-like pulses in the entire waveform and 
then looks for Sl-like pulses that precede the first S2-like 



pulse. Fig. [16| shows a typical events with SI and £2 and 
their hit patterns on the PMTs. 

The S2 pulse finding algorithm starts by applying a 
digital filter to the entire waveform to smooth out the 
high frequency components in order to facilitate the de- 
tection of the extent of pulses. It then searches the fil- 
tered waveform for regions where the signal exceeds a 
threshold of 10 mV for at least 0.6 /is, a time interval 
large enough to contain at least one S2 pulse, and for 
which the preceding and following 0.2 /is have an average 
signal less than 5% of the maximum within the interval. 
Because of the long afterpulsing tails that follow large S2 
pulses the interval above threshold will often contain mul- 
tiple pulses. The algorithm then recursively searches for 
S2 like pulses within that interval. This is done by com- 
puting the extent of any potential pulse by starting from 
its maximum sample in the interval and going backward 
in the trace until either the signal drops below 0.1% of its 
maximum or the slope of the signal changes sign. This 
defines the left boundary of this pulse. The same pro- 
cedure is repeated going forward in the trace to find the 
right boundary. If the pulse found has a FWHM larger 
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FIG. 16: (Color online) A typical event showing SI and S2 
and their hit patterns on the two arrays of PMTs (from [41] ) 



than 0.4 jis it will be considered as an S2 pulse and its lo- 
cation and boundaries will be saved. The recursive search 
for S2 pulses then continues within the interval (exclud- 
ing the regions where any pulse might already have been 
found). The parameters of the four largest S2 pulses are 
kept. A key difference in the MATLAB-based analysis 
was the development of a custom double box filter [46] , 
with the width of the two boxes set to 10 and 100 digi- 
tizer samples, respectively. These choices correspond to 
the approximate width of SI and S2 signals. This fil- 
ter proved extremely adept at correctly identifying the 
smallest S2 pulses. 

The SI pulse finding algorithm searches the total wave- 
form for signal excursions of at least 12 mV above the 
baseline. The boundaries of the potential pulse are de- 
fined as the points where the signal drops below 0.1% of 
its maximum. If the 0.5 lis preceding the pulse and the 50 
ns following the pulse respectively have an average signal 
less than 0.5% and 4% of the maximum, and if the decay 
time of the pulse is faster than 200 ns, then its location 
and boundaries will be saved. The parameters of the two 
largest SI pulses identified are kept. 

Once all the pulses have been identified, the reduction 
enters its last stage, computing the quantities of inter- 
est. Each PMT waveform is integrated over the bound- 
aries of the SI and S2 pulses identified and the values are 
converted into photoelectrons using the PMT gains. For 
each pulse a number of quantities are computed: FWHM, 
height, mean arrival time, number of PMT coincidences, 
number of digitizer channels saturating, etc. The XY 
position of each S2 pulse is computed using two differ- 
ent position reconstruction techniques (see section VB). 
A few parameters of the entire waveform are also com- 



puted: its total area, the average and RMS of the baseline 
before any pulses, and finally the total SI and S2 signals. 



B. Event Position Reconstruction 

An important feature of the XENON10 detector is its 
ability to localize events in 3D, enabling background re- 
duction with fiducial volume cuts and multiple scatter 
events rejection. Due to the small spread of the drift- 
ing electrons by diffusion in LXe [47 , the proportional 
scintillation signal is localized to a region with the same 
XY coordinates of the interaction site. Near the edge 
of the TPC, field non-uniformities lead to systematic ra- 
dial displacements towards the center of about 1 mm. 
Charge focussing effects in the extraction region between 
gate mesh and anode lead to charge displacements of less 
than half the mesh pitch, i.e., < 1 mm. The event z- 
coordinate, along the electric field direction, is inferred 
from the time difference between the primary and sec- 
ondary signals, with a resolution which is < 1 mm given 
the precision of the time measurement. Events with en- 
ergy depositions at different z positions will produce more 
than one S2 signal and can thus be readily identified and 
rejected. 

To obtain the XY event localization, with a few 
mm resolution, much smaller than the dimension of the 
PMTs, the light pattern measured by the top PMT array 
must be de-convolved. Two light pattern recognition al- 
gorithms were developed and tested on XENON10 data: 
a minimum x 2 algorithm and a neural network (AW) 
based algorithm. 

The minimum x 2 algorithm compares the measured 
light pattern with simulated S2 signals, calculating the 
values 



M 



[Sk - g»(a:,t/)] s 



(2) 



for all simulated XY positions. In this formula Si and 
Si are the measured and simulated S2 signals (in num- 
bers of photoelectrons) on the i-th PMT, M is the total 
number of PMTs (48 in the top array), and g\ takes 
into account the uncertainties of both the measured and 
simulated signals. The minimum value of the x 2 ( x ^y) 
function is used as an estimator for the true position in 
the XY plane. If the statistical uncertainty in the sim- 
ulated S2 pattern is small, the main contribution to o\ 
is the fluctuation in the measured signal, which includes 
the statistical fluctuation on the photoelectric emission 
from the PMT photocathode a pe ^ and the fluctuation 
a 9i on the gain g it If the number of detected photoelec- 
trons is sufficiently large, a pe ^ can be approximated as 
<j pe ,i ~ \fS~i while a g . is measured for each PMT through 
its single photoelectron spectrum 



' gi 
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(3) 
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The second algorithm is based on a fast artificial neural 
network. Its adaptive properties enable estimation of the 
vertex position even with a degraded input pattern, e.g., 
due to non-functioning PMT channels. The speed of the 
algorithm is nearly independent of the dimension of the 
detector. 

Due to the cylindrical geometry of the active vol- 
ume, we chose an algorithm based on two different sub- 
networks, working simultaneously to estimate the polar 
coordinates, radius r and polar angle 0. The NNs are 
feed-forward multilayer perceptrons with two hidden lay- 
ers. We applied an hyperbolic tangent as activation func- 
tion for the neurons in the hidden layers, with linear out- 
put for the output neuron. The network was trained, 
using the backpropagation rule, on a set of 4 x 10 4 simu- 
lated events, with the same simulation code used for the 
X 2 algorithm. The events were randomly generated in 
the XY plane with a mean value of 2.5 x 10 3 pe/events. 
This number is sufficiently large to minimize statistical 
fluctuations in the total number of detected photoelec- 
trons, and is comparable to the typical S2 signals in the 
energy range of interest. 

The input for both NNs is a vector with the top PMT 
signals in number of photoelectrons, normalized by the 
sum of their signals. In this way the output of the NNs 
is almost independent of signal size. The MC-estimated 
error on the radial position r is ~ 0.7 mm for r < 80 mm, 
and rises to 1.2 mm near the edge of the active volume. 



Fig. 17 shows a map of the reconstructed interaction ver- 
tices produced by a Cs source placed near the cryostat. 



The source location is apparent. In fig. [18] the radial dis- 
tributions of the Cs-interaction vertices obtained though 
NNs and % 2 reconstruction algorithms are compared with 
the corresponding MC simulation. Clearly the NN dis- 
tribution is in better agreement with the MC simulation 
near the edge of the detector, where reflection of scintilla- 
tion light on the Teflon renders the vertex reconstruction 
more difficult. Both NN and y 2 distributions show very 
good agreement with the simulated data for r < 80 mm. 



C. Basic Quality Cuts 

A number of basic quality cuts were defined to re- 
move events that do not correspond to real energy de- 
positions in the detector, events with unphysical param- 
eters, events with misidentified features, etc.. Many of 
the basic quality cuts are based on the known physical 
properties of the direct scintillation light of LXe and on 
the characteristics of proportional scintillation of gaseous 
Xe. 

The basic quality cuts applied to the SI signal are on 
the number of PMT coincidences, on the width and on 
the mean arrival time. In order to clearly distinguish 
low energy SI signals from random single photoelectrons, 
the minimum coincidence requirement is set to 2 PMTs 
having a signal above 0.35 pe. The efficiency of this cut is 
shown in figure 19 In addition, the width at 50% height 
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FIG. 17: (Color online) XY map of the reconstructed vertices 
with an external 137 Cs source: the position of the source can 
be clearly inferred. 
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FIG. 18: (Color online) Data-MC comparison between the ra- 
dial distribution of the vertices in the case ofa 137 Cs source ir- 
radiation of XENON10: the MC distribution (greeen) is plot- 
ted with the real data distributions obtained through NN(red) 
and x 2 (plue). 



of the SI pulse is required to be larger than 20 ns and 
smaller than 150 ns (note that this is the width in the 
digitized trace, after the anti-aliasing filter of the FADC 
as described in section IV) and its mean arrival time is 



required to be shorter than 60 ns. These cuts ensure that 
the SI signal is not confused with noise or with small 
single electron S2 signals, for example. 

Multiple basic quality cuts are also applied to the S2 
signal. They consist of a cut to reject events with more 
than one S2 pulse, a cut on the width, a cut to remove 
events with ADC saturation and a cut on the proportion 
of the S2 signal seen by the top and bottom PMT ar- 
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rays. As WIMPs are not expected to scatter more than 
once in the detector, events with more than one S2 pulse 
are usually discarded. This is also necessary for other 
analyses, when looking at the spatial depende nce of the 
direct scintillation for example (section VI A3). In terms 
of pulse shape, the FWHM of the S2 signal is required to 
be larger than 0.5 /is but less than 1 /is, consistent with 
what is expected from the proportional scintillation gap. 
Events for which the S2 signal peaks at a higher volt- 
age than the input range of the FADC are also removed, 
although this effect is absent at energies lower than 50 
keVee. Monte Carlo simulations have shown that the ex- 
pected fraction of S2 light measured by the top (bottom) 
PMT array should be about 55% (45%) and is consistent 
with what is observed in calibration data. Consequently, 
events where the ratio of S2 light seen on the top and 
bottom arrays falls outside the expected range are also 
rejected. For example, very large SI signals from alpha 
particles misidentified as S2 pulses would be cut since the 
bottom array would measure a much larger proportion of 
the signal. 

Another cut, intented to remove events in which the 
XY position has been incorrectly inferred, is also applied. 
The cut requires the x 2 between the S2 signal distribution 
from data and the Monte Carlo simulated distribution at 
the position inferred from the Neural Network algorithm 
to be less than a threshold value. 

Finally, an additional important parameter, related to 
the quality check of the pulse finding algorithm, is defined 
as 



log 



51 + 52 
A-S1-S2 



(4) 



where A is the total area of the trace and SI and S2 are 
the total signals identified. This parameter is analogous 
to a signal to noise ratio since it will be very large when 
all the features of the waveform have been identified and 
are included in the total SI and S2 signals whereas it will 
be very low if some features have been not been identified. 



VI. THE XENON10 CALIBRATION DATA 
RESULTS 

A. Gamma Calibration 



To achieve a reliable and accurate detector calibra- 
tion, with minimum impact on WIMP search exposure 
time, the XENON 10 shield was designed to allow intro- 
duction of external calibration sources without exposing 
the detector cavity to outside air. The gamma calibra- 
tion sources used were 57 Co, 137 Cs, 60 Co and 228 Th. The 
gamma source calibration data are used not only to de- 
termine the energy scale based on SI (or S2) signals, but 
also to define the detector's response to electron recoils 
from background events. 



1. SI light yield and energy scale 

The direct scintillation (SI) light yield in XENON10 
depends on the event energy and position as well as on 
the electric field strength in liquid xenon. The operating 
field in XENON10 was 0.73 kV/cm. Measurements of 
the SI yield at this field were carried out with gamma 
sources. In the bulk region (r < 5 cm), the 662 keV 
photo-absorption peak of 137 Cs gives an average light 
yield of 2.2±0.1 p.e./keV. For 122 keV gammas from 57 Co 
source, the light yield is about 3.1 p.e./keV, for events 
with radial positions between 8 and 9 cm, as shown in 
The different light yield for different gamma 



Fig. 20 



ray energies is due to the different contribution of the 
electron- ion recombination [28]. We were also able to 
observe the 30 keV characteristic X-ray peak for events 
near the edge of the sensitive volume. 
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FIG. 19: (Color online) Probability of SI identification as a 
function of the number of required coincidences. 



FIG. 20: SI scintillation spectrum from Co gamma calibra- 
tion. 
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2. S2 yield and liquid xenon purity 

The proportional scintillation light S2 is proportional 
to the number of electrons liberated in the liquid by an 
ionizing event. The number of electrons which reach 
the liquid-gas interface depends strongly on the liquid 
xenon purity, which can be inferred from a measure- 
ment of the electron lifetime r [48 . The presence of 
impurities reduces the number of electrons produced in 
the liquid at time t=0 (N e (0)) according to the rela- 
tion: N e (t) = A^Oje - ^ 7 ". To determine the electron 
lifetime, we used the data from 137 Cs calibration, mea- 
suring the attenuation with drift time of the S2 signal 
associated with the full energy peak of 662 keV gamma 
rays. Throughout the dark matter search period an elec- 
tron lifetime longer than 2 ms was measured, correspond- 
ing to <C 1 ppb O2 equivalent impurity concentration in 
the LXe. A similar purity level was inferred during a 
subsequent calibration of XENON10 with gamma rays 
from neutron-activated Xe gas, carried out shortly after 
the first dark matter search run was completed. Fig. 21 
shows the electron lifetime measurement using S2 from 
the 164 keV gamma events uniformly distributed in the 
XENON10 active volume. The electron lifetime inferred 
from this calibration is 2.2 ± 0.3 ms, confirming the excel- 
lent purity level achieved in XENON10 during the dark 
matter search. 



p.e. on the top PMT array and 9.9 p.e. on the bottom 
PMT array respectively. 
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FIG. 22: (Color online) S2 from single electrons seen in the 
neutron calibration data, with the signal detected by the top 
PMT array on the x-axis, and bottom PMT array on the y- 
axis. These events were selected based only on pulse width 
(full width 0.2 fis - 3.0 /xs) plus a requirement that they would 
be at least 10 fis from the trigger region (either before or 
after) . 
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FIG. 21: (Color online) Determination of the electron life- 
time from calibration data. The band is the S2 signal mea- 
sured by the bottom PMT array as a function of drift time 
for fully- absorbed 164 keV gamma events from the activated 
Xe calibration of XENON10. 

As previously discussed, the amplitude of the S2 signal, 
for a given number of ionization electrons extracted from 
the liquid, depends on several factors: the gas pressure 
p, the gas gap x, the electric field E across the gap, the 
light collection efficiency of the PMTs, and the quantum 
efficiency Qe of the PMTs. The observed number of 
photoelectrons per electron drifting in the gas can be 
found from single electron pulses within the data. During 
the operation of the XENON10 detector, a class of pulses 
due to single electron emission from the liquid to the gas 
phase was observed (Fig. 22). From this observation, we 



3. Position- dependence of SI and S2 Signals 

The SI light collection efficiency (LCE) has a strong 
dependence on the event position, due to the effects of 
total internal reflection at the liquid-gas interface, the 
solid angle, the optical transmission of the grids and the 
teflon reflectivity. In order to obtain an accurate energy 
calibration, data from both Monte Carlo simulation and 
internal sources were used to study the position depen- 
dence of the signals. 

Fig. 23 shows the simulated SI LCE throughout the 



estimate that, on average, each electron produces 13.£ 



XENON10 detector, where we define LCE as the proba- 
bility of photons released in the detector to hit the pho- 
tocathode of a PMT; it does not include the QE of the 
PMTs. Some regions are light- but not charge-sensitive, 
most notably in the reverse field region between cath- 
ode and bottom PMT array. Some light sensitivity also 
exists in the xenon around and below the bottom PMT 
array, where stray light enters through openings between 
the Teflon cylinder and the PMTs. The primary light 
is predominantly detected by the bottom PMT array, 
due to total internal reflection at the liquid-gas interface, 
where the index of refraction changes from 1.61 [49] to 
1. The ratio of top to bottom SI signal ranges from 
~ 0.14 to ~ 0.3 from the bottom to the top of the drift 
region. Some light is absorbed by the meshes, by the 
PTFE walls, or by the liquid xenon. We calculated 92% 
transparency of the meshes and assumed 92% reflectivity 
for PTFE [28] . The absorption length of liquid xenon is 
taken as 100 cm and the scattering length as 30 cm [50] . 
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FIG. 23: (Color online) Simulated SI light collection efficiency 
throughout the XENON10 detector. 



FIG. 24: (Color online) SI light yield throughout the sensi- 
tive liquid xenon target from activated xenon calibration with 
164 keV gamma rays. 
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In addition to the simulations, we also used neutron ac- 
tivated xenon isotopes, primarily 131m Xe (164 keV gam- 
mas) and 129m Xe (236 keV gammas) to calibrate the de- 
tector's SI response throughout the entire active volume. 
The activated xenon isotopes were introduced into the 
XENON 10 detector following several days of neutron ac- 
tivation of natural xenon [FT]. Unlike external calibra- 
tion sources, activated Xe isotopes provided unique en- 
ergy and position calibration throughout the entire LXe 
sensitive volume. The excellent 3D position sensitivity of 
XENON10 allows event-by-event signal correction based 



on position. Fig 24 shows the SI light yield for 164 keV 
gamma rays throughout the sensitive volume. The data 
can be compared with the simulated SI response (Fig.[23|) 
to determine parameters such as PTFE reflectivity and 
liquid xenon absorption length, that is relevant for the 
light collection. The response of the S2 signal at dif- 
ferent XY positions was also studied using the 164 keV 
gamma calibration data (Fig. 25). The signal size dif- 
fers by about 50% from the edge to the center, due to 
the sagging of the meshes which changes the gap where 
proportional scintillation light is generated. 

We used the SI and S2 light response maps from ac- 
tivated Xe calibration as a function of event position to 
correct the signals in the final data analysis. Better en- 
ergy scale, electron/nuclear recoil discrimination and bet- 
ter energy resolutions were achieved after the position- 
dependent corrections of the two signals. 
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FIG. 25: (Color online) S2 light yield response at different 
XY positions from activated xenon calibration with 164 keV 
gamma rays. 



4- Combined energy scale 

A recombination-independent combined energy scale 
for electron recoils, as described in details in [52 , 53 , can 
be produced from the appropriate sum of light (SI) and 
charge (S2). This scale should be free of the nonlinearites 
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present in the Sl-based scale. The energy deposited by 
each event is determined by combining 51 and S2 signals 
as following, 



E = ( 



SI S2 



)-W u 



(5) 



where 51 and 52 are in units of number of photo- 
electrons (N pe ). a and j3 are experimentally determined 
parameters in units of N pe /photon and N pe /electron, re- 
spectively. Wtot is the average energy required to pro- 
duce either a scintillation photon or an ionization elec- 
tron in LXe. In XENON10, f3 is determined from the 52 
corresponding to the single-electron emission peak (see 
Fig. 22). The anti-correlation between 51 and 52 is due 
to electron-ion recombination fluctuation in LXe. Each 
recombined electron-ion pair will create one UV photon. 
Thus, a in equation [5] can be determined by (3 and the 
slope 6 in Fig. 
From the 164 



26 , according to the relation a = (3/ tan#. 



seV calibration peak, we obtain, 



W tot = 14.0 eV 



(6) 



in good agreement with a study in a small LXe de- 
tector [54]. The energy resolution of the XENON10 de- 
tector was investigated with gamma ray sources ( 5T Co, 
22 Na, 13T Cs, 228 Th) covering the energy range between 
122 keV to 2.6 MeV. An example, from 137 Cs 662 keV 
gamma rays, is shown in Fig [27] For comparison, we 
also plot the energy resolution, in Fig. [28] obtained by 
using only 51, only 52 and the sum of these two signals. 
At 1 MeV, the resolution from the combined energy mea- 
surement is about a factor of seven better than that from 
51 alone, and a factor of three better than that from 52 
alone. The energy determined from the combination of 
SI and S2 signals show a much more linear response than 
that based on SI or S2 alone (see Fig. 29). 



B. Neutron Calibration 

To understand the XENON 10 response to nuclear 
recoils, a neutron calibration was performed using 
a 3.7 MBq 241 AmBe source, emitting ~220 neu- 
trons/second. The calibration was done by exposing the 
XENON10 detector to the source for approximately 12 
hours, with a live time fraction of 0.92. The 241 AmBe 
source (attached to a steel rod) was inserted through a 
7 mm diameter hole in the XENON10 shield. The source 
was positioned next to the cryostat, between two 5 cm 
thick Pb bricks used to block high energy (a few MeV) 
gamma rays produced by the source. 

The emitted neutrons have energies ranging from 
0.1 MeV to 11 MeV, with a mean at 4.3 MeV. The cal- 
ibration data were recorded at a constant rate of 6.5 Hz 
during the exposure. The trigger setup was the same as 
used for the WIMP search, with the addition of a high en- 
ergy veto to reject events with energies above ~120 keV. 
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FIG. 26: (Color online) Anticorrelation between SI and S2 
for 164 keV and 236 keV 7 rays from activated xenon isotopes 
( 131m Xeand 129m Xe). 
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FIG. 27: Combined energy spectrum for single-scatter events 
from 662 keV gamma rays interacting within the 5.4-kg fidu- 
cial mass of XENON10. Insets: SI and S2 distributions for 
662 keV photo-absorped events (events within 2-a around the 
662 keV peak in the combined energy are selected). The en- 
ergy resolution at 662 keV is 10.7%, 5.4% and 2.5% for SI, 
S2 and the combined energy, respectively. 



Figure 30 shows the log 10 (S2/S1) vs. SI distribution 



from the neutron calibration run, after applying the qual- 



ity cuts discussed in Section \VC\ Two regions are clearly 
distinguished in Figure 30 region a) which defines the 
nuclear recoil band corresponding to single elastic scat- 
ters, and region b) which corresponds to inelastic neu- 
tron scatters with 129 Xe, which produce 40 keV gamma 
rays. Neutron inelastic scatters with 131 Xe will produce 
80 keV gamma rays, however the data do not show a 
peak around this value because of the events rejection by 
the high energy veto. The nuclear recoil band is used to 
determine XENON 10 discrimination power, described in 
Section IVTCl 



17 




100 200 300 400 500 600 700 800 900 1000 

Energy [keV] 




°0 



20 



40 



60 



80 



100 



120 



S1 [keVee] (2.2 p.e./keVee) 



FIG. 28: (Color online) Energy dependence of resolutions in 
liquid xenon for SI (scintillation) only, S2 (ionization) only 
and combined (scintillation plus ionization) energy spectra, 
obtained for gamma rays below 1 MeV in XENON10. 
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FIG. 29: (Color online) Number of scintillation photons (N p h) 
and ionization electrons (N e ) per unit energy from gamma 
rays interaction in liquid xenon via photo-absorption for dif- 
ferent stopping power (dE/dx) at the operating field (0.73 
kV/cm) of XENON10. Values of dE/dx corresponding to dif- 
ferent energy of electrons (from photo-absorption of gamma 
rays) can be obtained from the ESTAR database [55]. The to- 
tal number of quanta (photons plus electrons) per unit energy 
does not depend on the energy of gamma rays. 



C. Nuclear/Electron Recoil Discrimination 

The success of XENON10 as a dark matter detec- 
tor hinges in large part upon its ability to discriminate 
electronic recoils from nuclear recoils, which in turn re- 
quires adequate definition of the detector response to 
such events based on calibration data. In addition to the 
formation of excitons, recoiling particles will produce a 
population of ionized electrons, many of which promptly 
recombine with their parent ions. Under an applied elec- 



FIG. 30: \og 10 (S2/ SI) vs. SI distribution from neutron in- 
teractions in the XENON10 detector. The data show two 
different regions: region a) which corresponds to recoils from 
elastic scatters and defines the nuclear recoil band, and region 
b) which corresponds to recoils from inelastic scatters. 



trie field, the relative number of recombining electrons 
decreases. However, because nuclear recoils have a char- 
acteristically higher ionization density than electronic re- 
coils, fewer electrons escape recombination from recoiling 
nuclei than electrons, for a given energy and drift field. 

Electrons which recombine contribute to the prompt 
scintillation signal (SI), while those which escape recom- 
bination are drifted to the anode in the gas and produce 
the proportional signal (S2). The relative strength of re- 
combination for a given event can be measured by the 
ratio S2/S1, and hence this parameter can be used to 



discriminate between recoiling species. Figure [31] shows 
the behavior of Logio(S2/Sl) as a function of energy , 
in unit of keV electron-equivalent (keVee), for popula- 
tions of both recoil types, called the electronic and nu- 
clear recoil bands, or ER and NR bands, respectively. 
The main purpose of such ER and NR calibrations is to 
identify a region in Logio(S2/Sl) vs. SI space, called 
the WIMP acceptance window, which should be nearly 
free of ER events while covering a significant portion of 
the NR band. The lower bound of this window along 
the horizontal axis is determined by the detector's SI 
threshold, and the corresponding upper bound is chosen 
to maximize the potential integrated WIMP rate while 
minimizing the effects of anomalous backgroun d event s 
which occur mostly at higher energies (see section VII B ) . 



The choice of bounds along the vertical axis are discussed 
here. 

Data were taken with the 137 Cs source throughout 
all of November 2006, and intermittently from Decem- 
ber 1 through February 14, 2007, accumulating a to- 
tal of ~ 2100 events (after quality and fiducial cuts) in 
the WIMP acceptance energy window, 4.4 p.e. < SI 
< 26.4 p.e. Fluctuations in Logio(S2/Sl) over most of 
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FIG. 31: (Color online) The electronic and nuclear recoil 
bands shown in Logio(S2/Sl) vs. SI space. 



this range are dominated by recombination fluctuations, 
until the lowest energies where uncorrelated statistical 
fluctuations take over. The width of the electronic re- 
coil band is very important for gamma rejection because 
the two bands overlap. Due mainly to the non-uniform 
SI response at different locations in the fiducial region, 
applying spatially-dependent corre ctions t o SI based on 
the 131m Xe calibration (see section VI A 3 ) improves the 



overall SI resolution and thus helps to reduce the vari- 
ance of the bands. Data with the AmBe neutron source 
were taken on December 1, 2006 for approximately 12 
hours, accumulating a total of about 260,000 events. The 
energy dependence of both bands makes it difficult to pre- 
cisely measure the discrimination power in the absence of 
extraordinarily large calibration datasets. In an effort to 
remove this energy-dependence, a one-dimensional trans- 
formation that "flattens" the ER band is applied to all 
data. The ER band is broken up into 1 keVee-wide, ver- 
tical slices in SI. For each, a Gauss fit is applied to the 
Logio(S2/Sl) spectrum. The mean of each fit now rep- 
resents the center of the ER band in that particular bin. 
A high-order polynomial is fit to the Gauss means, which 
provides an analytic form for the ER band centroid as a 
function of SI, and is subtracted from every data point in 
both bands. This procedure flattens the ER band (and 
to a large extent, the NR band as well), and introduces 
a new parameter, ALogio(S2/Sl), which represents the 
distance from the ER centroid in Logio(S2/Sl) space. 
Figure 32 shows the bands in ALogio(S2/Sl) space. 

Although the energy dependence of the ER band cen- 
troid has been removed, the NR band centroid and width 
still change with energy. Again, the flattened bands are 
broken up into vertical SI slices, only this time more 
coarse binning is used — seven bins in the WIMP en- 
ergy region of interest (ROI) — in order to maximize the 
statistics in each slice, and a Gauss fit is applied to the 
ALogio(S2/Sl) spectrum of both bands. One such slice is 
shown in Figure [33} for the range 13.4-17.2 keVr (1 keVr 
= 1.1 pe, according to [11 ) of nuclear recoil energy. The 
WIMP acceptance window is defined to lie in the range 
(/i — 3<j) < ALogio(S2/Sl) < /i, where \i and a are the 
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FIG. 32: (Color online)The bands in figure pi] have been 
transformed to show the distance in Logio(S2/Sl) space from 
the ER band center, giving the new discrimination parameter, 
ALogi (S2/Sl). The vertical lines indicate the WIMP region 
of interest (ROI). 



mean and sigma from the NR band Gauss fits, respec- 
tively. The Gauss fits were performed only to define the 
window bounds; the NR acceptance, A nri was calculated 
by counting the number of AmBe events that fall within 
this window, for each energy bin. 
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FIG. 33: (Color online) Distributions of ALogi (S2/Sl) for 
nuclear and electronic recoils in the range 13.4-17.2 keVr. The 
WIMP acceptance window in this particular energy range is 
denned by the blue, shaded rectangle which is between fi and 
\i — 3a of the NR band. 



The shape of the ALogio(S2/Sl) fluctuations in the 
ER band are "empirically" Gaussian; that is, with 
the statistics available, they appear consistent with 
a Gaussian distribution. As previously stated, the 
ALogio(S2/Sl) spectrum is dominated by recombina- 
tion fluctuations, which are poorly understood, and thus 
more cannot be said in the absence of a larger calibra- 
tion dataset. We calculate the predicted ER rejection in 
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TABLE II: The nuclear recoil acceptance A nr , the efficiency e cu t 
of remaining nuclear recoil events after the GammaX cuts (see 
sect ion [Vll| , and the electron recoil rejection efficiency R er for each 
of the seven energy bins. The predicted number of leakage events, 
Nieak, is based on R er and the number of background events, N ev t, 
in each energy bin, for the 58.6 live-days WIMP-search data. Errors 
are the statistical uncertainty from the Gaussian fits on the electron 
recoil ALogio(S2/Sl) distribution. 

(keV) A 1 - R er N evt Ni eak 

(io- 3 ) 

i+U.7 



4.5 - 6.7 
6.7- 9.0 
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the case that ALogio(S2/Sl) fluctuations are Gaussian. 
That is, we use the Gauss fits to the ALogio(S2/Sl) spec- 
trum in each of the seven energy bins to determine the 
energy-dependent discrimination power. The results are 
shown in Table [II] and Figure [34] Additionally, the ex- 
pected number of background events in the WIMP accep- 
tance window, Ni ea k, is shown, calculated on the basis of 
predicted rejection and ba ckgro und rate in the 58.6 live- 
days exposure (see Section VII). 



Figure [35] shows a decomposition of the ER band vari- 
ance, using assumed statistical and instrumental fluctu- 
ations. It is quite evident that uncorrelated statistical 
and instrumental fluctuations cannot alone account for 
the observed degree of variance. Unfortunately, a model 
does not yet exist that successfully predicts recombina- 
tion fluctuations in noble liquids, and hence more mea- 
surements are needed to reach a better understanding of 
the subject. 
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FIG. 35: (Color online) Decomposition of the ER band vari- 
ance. The recombination fluctuations are inferred by compar- 
ing the expected statistical and instrumental fluctuations to 
the full observed band variance. 
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FIG. 34: The ER rejection as a function of SI for 
ALogio(S2/Sl) < \i. The rejection improves at lower en- 
ergies, to better than 99.9% in the range 2-3 keVee. 



The observed trend of the ER rejection power with en- 
ergy is unexpected. If recombination fluctuations were 
flat at all energies, or if the band widths were dominated 
by binomial fluctuations from light collection, photoelec- 
tron emission, etc., one would expect the band widths to 
grow at low energies, and hence the ER rejection power 
would deteriorate. The opposite is observed, and is due 
to two factors. First, the ER and NR bands themselves 
diverge slightly at lower energies, as first measured in [4]. 
Second, the width of the ER band does not grow at 
lower energies but instead remains relatively constant. 



VII. THE XENON10 DARK MATTER SEARCH 
DATA 

The XENON10 Dark Matter Search data covered the 
period between October 6th 2006 and February 14th 



2007. Figure 36 shows the accumulation of measurement 
livetime, interrupted by calibration runs. The global 
trigger rate was about 2.6 Hz with more than 92% live- 
time. A total of 74.7 live-days of data were collected. Of 
this, 16.3 were selected to be analyzed in order to de- 
fine the event selections and cuts. The remaining 58.4 
live-days were embargoed in order to carry out a blind 
analysis. The XENON10 results from this blind analysis 
have been published in [TTJ [12], and are referred to as 
blinded WIMP search data in the following. 



Background reduction with fiducial volume 
selection 



Background events in XENON 10 are dominated by 
electron recoils from radioactivity of detector materials. 
Due to the self-shielding of LXe, the background event 
rate falls dramatically from the edges to the center of the 
target. Based on Monte Carlo simulations and unblinded 
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FIG. 36: (Color online) Running days vs live-days from 
XENON10 data-taking. 



background data, a central fiducial volume was chosen to 
minimize background while maintaining the largest pos- 
sible target mass. The fiducial volume was pre-defined 
for the blinded analysis with a radius of 80 mm and a 
drift time between 15 and 65 /is, corresponding to a fidu- 
cial mass of 5.4 kg. Figure [37| shows the fiducial volume 
overlaid on the background rate distribution of single re- 
coil events between 2-30 keVee (2.2 pe/keVee). Within 
the chosen fiducial volume, the overall background rate 
is reduced by almost a factor of 10, as shown in Figure 

M 

A comparison between the measured background rate 
and the projected background from detailed Monte 
Carlo simulations, shown in Figure [39| reveals a dom- 
inant contribution from the stainless steel vessels, fol- 
lowed by approximately equal contributions from PMTs, 
feedthroughs, and Teflon. The contribution from intrin- 
sic contamination of 85 Kr, Rn, 136 Xe in the liquid xenon, 
or nuclear recoils from (a,n) reaction or fission of 238 U 
in the surrounding environment is subdominant or neg- 
ligible in the entire WIMP search run. Details of the 
simulations and understanding of background contribu- 
tions are discussed in Appendix [Aj 



B. Definition of anomalous events and cuts 

During the analysis of both gamma calibration data 
and unblinded background data from XENON10, a class 
of anomalous single scatter events with smaller S2/S1 
ratio than that of normal single scatter electron recoil 
events were observed. These events produce a non- 
Gaussian tail extending into the acceptance region for 
nuclear recoils and to even lower 52/51 ratios. Most 
of these anomalous events are likely due to the coinci- 
dence of two gamma-ray interactions, where a single in- 
teraction occurs in the sensitive volume and at least one 
other interaction occurs in a region of the detector that 
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FIG. 37: (Color online) 2-30 keVee (based on SI and 2.2 
pe/keVee) event rate distribution in the fiducial volume as a 
function of the drift time and radius for the WIMP search 
run. The fiducial volume cut used in the final analysis is 
pre-defined based on the gamma calibration and unblinded 
WIMP search data and is indicated with the dashed lines. 



3L 30 




No Fiducial cut 

R<80 [mm] and 1 0<dt<65 [\i s] 




1000 



1500 



2000 2500 3000 
Energy [keVJ 



FIG. 38: (Color online) Overall XENON10 background en- 
ergy spectrum before (black) and after (red) the fiducial vol- 
ume cut. Here we use the combined energy scale. 



is light-sensitive but not charge-sensitive. Such events, 
dubbed Gamma-X, may comprise two Compton scatters 
of a high-energy gamma-ray, or, e.g., coincident gamma- 
rays from a cascade, such as 60 Co from the stainless steel 
cryostat. There are a few regions of the XENON 10 detec- 
tor that are light-sensitive but not charge-sensitive, such 
as the region between cathode and bottom PMTs or the 
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FIG. 39: (Color online) The measured XENON10 background 
event rate (black) compared to Monte Carlo simulations (red) 
in the 5. 4- kg fiducial mass. The energy is based on com- 
bined energy scale. The contribution from different compo- 
nents from the main detector and shield are shown in different 
colors. FTs represent feedthroughs. Poly represents polyethy- 
lene in the shield. IC and OC represent internal cryostat and 
outer cryostat separately. 



region around and below the bottom PMT array, as was 



shown in Fig. 23 



For illustration, Fig [40] shows the ALogi (S2/Sl) vs. 
energy event distribution of the blinded WIMP search 
data. Prior to a final set of cuts, designed to remove these 
anomalous events, a total of 22 events is observed in the 
WIMP acceptance region. Only some of these events are 
expected from the statistical leakage of the electron recoil 
band. The spatial distribution of all events in the energy 



region of interest is shown in Fig 41 

Two cuts, based on the SI hit pattern, were identified 
to remove these anomalous events. One cut is based on 
the asymmetry of SI between the top and bottom array 
of PMTs, i.e. p asy = (Sl t - Sl b )/(Sl t + Sl b ). p aS y 
is larger for anomalous events with some scatters in the 
region close to the top PMTs, but not in the sensitive 
LXe target. This cut also removes most of the events 
with misidentified SI in software, e.g. a small S2 signal, 



or a micro-discharge pulse in the gas. Figure [42] shows 
the distribution of this parameter versus energy. 

A second set of cuts is based on the SI signal distri- 
bution on the top and bottom PMT arrays. We define 

a parameter Prms = \J\ XX^l* — ^l) 2 fr° m the Root- 
Mean-Square (RMS) value of the 5 (for top) or 10 (for 
bottom) PMTs receiving the largest hits. Anomalous 
events tend to give non-uniform Sl-signal distributions, 
e.g., if one of the interactions occurs near the bottom 
PMT array. Figure [43] shows the effect of this cut for the 
bottom PMTs. 

The above software cuts also remove a small fraction of 
genuine nuclear recoil events. By comparing the AmBe 
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Energy [keVee] (based on 2.2 phe/keVee) 

FIG. 40: (Color online) Event distribution in the 
ALogi (S2/Sl) vs. Energy space for the 58.4-day WIMP 
search data. The crosses indicate the events that "leak" into 
the nuclear recoil acceptance region, which is within the mean 
and 3cr of the nuclear recoil band. Red circles indicate the re- 
maining events after all software cuts, discussed in the text. 
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FIG. 41: (Color online) Event distribution in the 58.4-day 
WIMP search data. The crosses indicate the events that 
"le ak" into the nuclear recoil acceptance region, as shown in 
Fig [40] in the energy window (2-12 keVee) of interest. Red 
circles indicate the remaining events after all software cuts, 
discussed in the text. 



nuclear recoil data before and after the software cuts, the 
estimated cut efficiency e cut was calculated and shown in 
Table HI 



Appendix A: The XENON10 Backgrounds 

The particles potentially contributing to the back- 
ground of the XENON10 detector are alphas, betas, gam- 
mas, and neutrons. Neutrons scattering with Xe nuclei 
yield nuclear recoils, possibly mimicking a WIMP sig- 
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FIG. 42: (Color online) A cut based on the asymmetry of the 
SI signal, p asy , as discussed in the text. The crosses indicate 
the events that "leak" into the nuclear recoil acceptance re- 
gion, as shown in Fig |40| The two vertical lines are the energy 
window of WIMP search. Red circles indicate the remaining 
events after all software cuts. The solid black curve is pre- 
defined to keep 98% electron recoil events. It removes 3 out 
of 22 "leakage events" in the energy window of interest (2-12 
keVee) . 
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FIG. 43: (Color online) Cuts based on pRMSb for the bottom 
PMTs, as defined in the text. The crosses indicate the events 
that "leak" into the nuclear recoil acceptance region, as shown 
in Fig |40| The two vertical lines are the energy window of 
WIMP search. Red circles indicate the remaining events after 
all software cuts. The solid black curve is pre-defined, which 
varies from 99% at 2 keVee and 90% at 20 keVee cut efficiency 
and fitted by the third-order polynomial function. It removes 
9 out of 22 "leakage events" in the energy window of interest 
(2-12keVee). 



nal. Unless they are detected as multiple interactions 
in the sensitive volume, neutrons constitute an unre- 
jectable background in the nuclear recoil band. Gamma 
and beta rays mostly populate the electron recoil band 
and constitute the most prominent background of the 
experiment. Alphas can contribute to background in the 
WIMP search region as they give rise to neutrons via a-n 



reactions in the detector materials. 



Gamma and Neutron Backgrounds from Data 
and MC Simulations 



Gamma and neutron backgrounds from the environ- 
ment were studied to optimize the shield, and have been 
discussed in section III After shielding, the dominant 



backgrounds in XENON10 result from radioactive con- 
taminants in the detector and shield materials. These in- 
ternal backgrounds were studied with material screening, 
using gamma-ray spectroscopy in a low count rate facil- 
ity, with Monte Carlo simulations, and with the analysis 
of delayed coincidences in decay chains of radio-nuclides 
inside LXe. Simulations were based on the GEANT4 
simulation toolkit [56 , where the complete geometry of 
the detector and surrounding material was included. 

As detailed in Table |TJ the shielding material and many 
detector construction components were screened with 
HPGe gamma spectrometers at the LNGS and SOLO 
screening facility [57]. The radioactive contamination of 
the screened material are lis ted in Table [I] (for the shield) 



and in the first row of table III (for the PMTs + bases). 



However, since XENON 10 was a prototype experiment, 
the radioactivity of some detector components was not 
known. We used the minimum y 2 method to extract the 
individual activities, by comparing the measured energy 
spectrum with the sum of the energy spectra for each 
radioactive decay chain within each material: 



N 



di — Y^7=i Fij ( u j Uj •> thj Thj , kj Kj , coj C oj , csj C sj 



where n and N are the first and last bin of the energy 
spectra, defining the fit range. Here we use a bin width of 



3.75 keV, 



50, and N = 450. dj is the rate in the 



Ath 



bin of the measured energy spectrum, and is the cor- 
responding expected rate from simulated backgrounds in 
material j; m = 41 is the total number of simulated ma- 
terials; Uj, thj, kj, coj and csj are the scaling parameters 
related to nominal activities of material j in 238 U (Uj), 
232 Th (Thj), 40 K (Kj), 60 Co (Coj), and 137 Cs (C Sj ), 
respectively. The minimization is done by the MINUIT 
routines in ROOT 05]. 

The results of this fit procedure are summarized in 



table |III| The main component to the background energy 
spectrum of XENON10 originates from a 137 Cs peak at 



662 keV. Fig. 44 shows the 662 keV peak decreasing at 
smaller radii, a clear indication that it is not diluted in 
the LXe. A study of the spacial distribution of events 
in this peak shows that the 137 Cs source is all around 
the detector edges. Thus a surface contamination of the 
anthropogenic 137 Cs has to be assumed, either on the 
inner cryostat can or on the Teflon. The best fit results 
localize the 137 Cs on the Teflon, with a total activity of 
5.73±0.57 Bq. 
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TABLE III: Radioactive contaminations of material used in the construction of the XENON10 detector, as obtained from 
screening (PMTs, bases and polyethylene shield) and by comparing the low-background data to MC simulations. 
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0.23 
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37.4 kg 
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470 
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FIG. 44: (Color online) Energy spectrum of background 
events in XENON10 with different radial cuts, with a focus 
on the 137 Cs peak. 



The main source of neutrons in XENON10 are (a,n) 
and spontaneous fission reactions from 238 U and 232 Th 
in the detector and shield material. Another source of 
neutrons are spallation and photo-nuclear reactions of 
cosmic ray muons in the rock and shield. As explained 
in Section |III[ the neutrons from outside the shield are 
stopped or thermalized by the polyethylene. 

Using the 238 U and 232 Th activities as given in Ta- 
ble [IllJ the energy spectra and number of expected neu- 
trons from (a,n) and spontaneous fission reactions in 
each material has been calculated with the modified 
SOURCES4A code [58]. These are given in Table [IV 
(second column). 

The neutrons are then propagated into the sensitive 
region using the GEANT4 code and detector geometry. 



In figure [45] the predicted single scattering nuclear recoil 
energy spectrum from the neutrons from all the mate- 
rial considered in the simulation is shown, together with 
the individual contributions from all the material. The 
number of predicted single nuclear recoils in the WIMP 
search region is shown in Table IV (third column), along 
with the total number of single nuclear recoils (fourth 
column) . 



The total rate of single nuclear recoils in the energy 
region from 2 to 12 keVee is expected to be 1.67xl0 -3 
event/kg/day, i.e. for an exposure of 5.4kgx58.6 days, 
the total number of nuclear recoils expected from the 
contamination of the XENON 10 detector material is less 
than 0.53 events [TT] . 
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FIG. 45: (Color online) Monte Carlo simulation of total nu- 
clear recoil spectrum expected from the neutrons produced in 
all the material considered (black) , along with the individual 
contributions from each material. 



2. Intrinsic backgrounds in LXe 

Delayed coincidence analyses to determine the Radon, 
Thoron, U/Th and 85 Kr concentration in the LXe have 
been carried out by looking at the specific decay signa- 
tures. For Radon, the consecutive beta-alpha d ecay s oc- 



curring in the decay of 214 Bi (shown in Eq. |A2| are 
required to occur within the same waveform, in the same 
time order as in the decay scheme and with the proper en- 
ergies (E/3 < 3000 keV and E a > 6500 keV). With these 
cuts, the residual accidental coincidences are negligible 
and the events surviving the cuts are taken as real 214 Bi 
decays. The resulting Rn concentration in LXe, consid- 
ering the detection efficiency and the efficiency of above 
cuts, is (59±2)/iBq/kg. For 220 Rn arising from the 232 Th 
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TABLE IV: Neutron flux es (second column) calculated with the modified SOURCES4A code, using the activities in 238 U and 



Th given in table III Column 3 lists the number of single nuclear recoils expected in the WIMP search region of interest, 
inside the fiducial volume, for each material. Column 4 lists the total number of single nuclear recoils expected in the whole 
volume. 



Material 


Flux 




nuclear recoil total nuclear recoil 




[neutron/sec] 


evt/kg/day 


evt/kg/day 


Teflon 


9.60x10" 


-6 


5.34xl0 -4 


3.69xl0" 4 


Poly 


5.22x10" 


-6 


2.38xl0 -6 


7.32xl0~ 7 


OC 


5.64x10" 


-5 


2.82xl0 -4 


1.88xl0 -4 


IC 


4.38x10" 


-5 


5.57xl0" 4 


4,02xl0~ 4 


PMT 


2,41x10" 


-8 


1.13xl0 -4 


5.16xl0 -4 


Base 


3.81x10" 


-6 


1.80xl0 -4 


6.39xl0" 5 



chain, the time window for the events is much shorter, 
namely 3 /is versus the maximum time acquisition win- 
dow, i.e. 80 /x, and the energy of the two particles has to 
be Ep < 2400 keV and E 1 >7000 keV. The resulting con- 
centration of 220 Rn in LXe is (4.7±0.2)/iBq/kg. For both 
Radon and Thoron the events are distributed mainly in 
the outermost LXe region, suggesting that both radioiso- 
topes emanate mainly from the surface of the teflon vessel 
or inner stainless steel chamber. 

The delayed coincidence technique has also been used 
to estimate the Krypton concentration in LXe, by tagging 
the beta-gamma decay mode (0.434 % branching ratio) 
of the 85 Kr diluted in the liquid. 85 Kr, which has a con- 
centration about 10 -11 in Kr, is a beta emitter that pro- 
duces background in the LXe volume. XENON 10 exper- 
iment requires the Kr/Xe contamination to be less than 
10 ppb, which gives about 0.2 dru electron recoil back- 
ground. To have this low Kr/Xe contamination level, the 
Xe gas used for XENON10 was purified by the Spectra 
Gases company [59] in several passages through a large 
cryogenic distillation column. Figure 46 shows the en- 
ergy spectra of the gamma and beta particles selected 
by the cuts. A total of 59 delayed-coincidence events 
were found in the entire WIMP-search data. Taking into 
account the detection efficiency of these events, we de- 
termine an 85 Kr activity of (1.0±0.06) mBq/kg, which 
translates into (5.0±0.3)ppb of natural Kr/Xe contam- 
ination for typical atmospheric 85 Kr abundances. This 
value is compatible with the value quoted by the Xe gas 
supplier [59] and is sufficiently low for XENON10. 
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FIG. 46: (Color online) Gamma and beta energy spectra of 
the 85 Kr - 7 decay (0.434% BR) from the XENON10 data, 
after delayed coincidence cuts. 
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